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Abstract 


We  have  made  a  quantmn  mechanical  perturbed  stationary  state 
calculation  of  the  rate  of  recombination  of  oxygen  atoms  in  the  presence 

of  an  oxygen  molecule  as  third  body.  The  recombination  is  considered  to 

3  ^ —  * 
take  place  into  highly  excited  vibrational  levels  (v  =  9,10)  of  the  A  2_ 

state  of  0  .  There  are  no  adjustable  parameters  in  the  calculation:  for 

the  interaction  potential  we  have  used  the  results  (for  0-  X  Yl  p  ^   obtained 

in  a  previous  calcvilation,  which  were  obtained  by  comparison  with  experimental 

measurements  of  vibrational  relaxation  times  in  oxygen.  At  temperatures  of 

the  order  300  -  3000°K,  the  calculated  recombination  coefficient  into  these 

two  vibrational  levels  is  of  the  order  2  -  10  x  10"-^-^  cm  /sec,  which  agrees 

to  within  better  than  a  factor  of  3  with  recent  experimental  measurements. 
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1.  Introduction 

a.  Statement  of  the  Problem 

Two  atoms  can  recombine  to  form  a  molecule  in  the  presence  of  some 
third  body  such  as  an  atom  or  molecule  (or  even  a  wall,  but  we  shall  not 
consider  this  case),  or  else  radiatively.  In  practice,  three-body  non- 
radiative  recombination  is  an  important  and  interesting  process  which  has 
considerable  practical  importance  in  the  ionosphere,  in  shock  wave  relaxation, 
and  in  combustion;  in  fact  three-body  recombination  is  the  prototype  of  an 
important  class  of  chemical  reactions. 

For  definiteness,  we  shall  examine  the  recombination  of  two  oxygen 


atoms  in  presence  of  an  oxygen  molecule  - 

(1.1)  0   (^P)    +  0   (^P)    +  0      (X^  H'J   ^    0     (A^  Yl*^    +0      (X^  XI  "). 
which  we  shall  write  as 

(1.2)  A  +  B  +  C  ->  Ab'    +  c"     . 

The  primes  indicate  that  the  molecule  AB  is  in  general  formed  in  an  internally 
excited  state,  and  the  third  body  C  may  also  carry  away  some  internal  excitation 
after  impact. 

It  must  be  stressed  that  the  present  calculation  is  quite  schematic; 
thus   the  detailed  labeling  Op (A  ^  )  is  not  to  be  taken  too  seriously.  In 
fact,  we  have  chosen  this  state  because  the  vibrational  levels  of  this  state 
are  known  right  up  to  the  dissociation  limit'-  -" .  The  purpose  of  the  present 
calculation  is  to  sketch  the  general  order  of  magnitude  of  a  three-boc^ 
recombination  process,  rather  than  to  make  detailed  calculations  for  one 
specific  process. 

We  con»ider  the  case  of  a  molecule  AB  formed  in  a  highly  vibrationally 
excited  state,  but  do  not  consider  the  third  body  C  to  have  any  internal  ex- 
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citation  either  before  or  after  impact.  On  theoretical  grounds  ^  '  -'  ve 
believe  that  the  raolecvile  AB  is  very  likely  to  be  formed  in  a  highly  excited 
vibrational  state,  because  in  this  way  a  minimum  of  energy  has  to  be  trans- 
ferred into  kinetic  energy  of  C  from  the  'continuum'  (i.e.,  imbound)  vibration 
of  A  and  B,  and  it  is  well  known  that  the  coupling  of  translational  and  vibrational 
energy  is  weak,  so  that  energy  transfer  between  these  two  modes  of  motion  is 
difficult'-  -' .  It  should  however  be  noted  that  thus  far  there  appears  to  be  no 
experimental  evidence  for  the   existence  of  such  highly  excited  vibrational 
states'-  -'.  The  observations  would  be  very  difficiilt,  but  one  should  make 
fiu:'ther  attempts  to  make  the  measurements,  as  well  as  to  calculate  the  rate  of 
degradation  of  vibrational  energy  in  collisions. 

Let  us  examine  the  implications  of  this  model  somewhat  further.  We 
may  discount  electronic  excitation  of  either  AB  or  C  from  adiabatic  considerations 
(except  in  somewhat  special  situations  that  do  not  concern  iis  here).  Also, 
except  in  the  special  case  of  exact  resonance,  we  may  disregard  any  vibrational 
excitation  of  C  after  the  collision,  on  account  of  the  well-known  difficulty 
of  transferring  relatively  large  amounts  of  energy  between  translation  and 
vibration  L-" .  There  remains  the  problem  of  rotational  motion.  It  is  clear 
that  in  most  actual  collisions  there  will  be  relative  angular  motion  of  AB  and 
C,  and  the  effects  of  this  should  certainly  be  examined.  However,  it  is  well 
known'-  -^  that  transfer  of  energy  between  translation  and  rotation  is  easy,  and 
thus  we  feel  that  we  can  obtain  an  initial  physical  understanding  of  the  situation 
by  merely  examining  a  collinear  collision,  in  which  rotational  motion  is  dis- 
regarded. For  definiteness  we  examine  recombination  into  the  two  highest 
vibrational  states,  v  =  9,10. 

b.  The  Present  Model 

We  consider  the  collinear  collision  of  two  atoms  A,  B  with  a  third 
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body  C.  The  free  atoms  A,B  are  regarded  as  in  the  vibrational  continuiora  of 
some  electronic  state  of  the  molecule  AB,  and  they  fall  into  a  boimd,  but  still 
highly  vibrationally  excited,  state  of  AB  (still  of  the  same  electronic  state 
of  the  molecule).  The  balance  of  the  'vibrational'  energy  of  A  and  B  is  trans- 
ferred into  kinetic  energy  of  the  third  body  C,  measured  relative  to  the  center 
of  mass  of  AB.  Thus  we  carry  out  a  calculation  of  the  transfer  of  energy  between 
vibrational  and  translational  degrees  of  freedom,  analogous  to  earlier  calculations 
of  vibrational  relaxation  times  in  gases '-'-'. 

This  is  an  adequate  model  within  its  very  obvious  physical  limitations 5 
but  in  practice  we  use  a  quantum  mechanical  perturbation  method  of  calculation, 
the  'method  of  pert\irbed  stationary  states',  which  is  adequate  for  vibrational 
relaxation  processes  (i.e.,  the  v  =  0  to  v  =  1  transition,  which  has  a  very 
small  transition  probability),  but  not  for  a  process  like  three-body  recombination, 
which  is  a  very  probable  process.  The  situation  is  the  following.  One  can 
calcvilate  reaction  rates  classically,  with  a  relatively  detailed  statistical 
model  but  only  a  crude  estimate  of  transition  probabilities'^'  *   ,  or  else  one 
can  carry  out  a  quantiim  mechanical  calculation  with  a  less  detailed  statistical 
model,  but  a  better  estimate  of  transition  probabilities'-  '  *  -• .  To  get  the 
complete  story,  one  needs  a  synthesis  of  these  two  approaches:  for  processes 
having  small  transition  probabilities  the  quantum  method  is  better,  while  the 
classical  approach  gives  more  information  for  processes  having  large  transition 
probabilities.  Nevertheless,  there  are  certain  features  that  cannot  be  treated 
adequately  by  a  simple  classical  model,  such  as  the  dependence  of  the  cross 
section  on  the  final  state  of  vibrational  excitation,  and  it  is  for  this  kind 
of  detailed  information  that  a  calculation  like  the  present  one  is  useftil. 


We  should  like  to  thank  Dr.  N.  Davidson  and  Dr.  J.C.  Keck  for  some  relevant 
remarks . 


Thus   it  is  necessary  to  examine  the  niimerical  resiilts  of  the  present  work 
very  carefully  as  a  basis  for  improving  one's  general  understanding  of  the 
problem. 

It  should  be  pointed  out  that  in  the  present  case  the  intermolecular 
potentials  are  known  even  less  well  than  for  the  calculation  of  vibrational 
relaxation  times'-  *  -"j  but  becai:ise  the  details  of  the  quantum  mechanical  cal- 
ciilation  are  obscured  by  the  large  transition  probability,  this  is  not  of  the 
same  importance  as  in  other  processes. 

2.  Calculation 

a.  Method  of  Perturbed  Stationary  States 

A  rather  thorough  discussion  of  this  method,  and  of  the  techniques  of 
solution  used  here,  has  been  given  previously'-  *  -^ .  Call  the  two  atoms  that 
will  combine  A  and  B,  and  the  third  body  C.  Then  we  use  the  following  coordinates, 
in  which  the  motion  of  the  center  of  mass  separates  off,  and  v/ith  which  there 
are  no  cross  terras  in  the  kinetic  energy  - 

r 


(2.1) 


p  =  Ia  -  Ib 


If  we  separate  off  the  motion  of  the  center  of  mass  in  the  coordinate  R,  the 
wave  equation  becomes 

(2.2)  (H  (p)  +  H  (r)  +  v'(p,r)  -  E)  I.(p,r)  -  Oj  j  -  f,b, 

1  —     e        -         J  — 

where 


(2.3) 


H^(p)  -  -  (tiV2M^)Vp  +  V°(p)  J  M^  =  M^  Mg/  (M^  +  Mg) 


P 
.2  /„..  ^^2 


H^(r)  -  -  (ii^/2M^)y^  +  V°(r)  :  M^  =  (M^+Mg)Mj,/(M^+Mg+Mp), 
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'f  and  'b'   stand  for   'free'  and   'bound'   states  respectively,  while    'i'  and 
'e'    stand  for   'internal'    (vibrational)  and   'external'    (translational)  motion. 

We  solve  the  wave  equation  (2,2)  by  the    'method  of  perturbed  stationary 
states',  in  terms   of  the  unperturbed   ffemiltonian  H     «=  H.(p)   +  H  (r)  - 


(2.1,) 


(H^  -  E)  Y.ip,T)   -  0  J  Y  j(p,r)  -  0.{p)  Y^(r) 

00 

Y   .(r)    'JZ         ^l/r)  X(V^   P<7    ^'=°^  ^i^> 
'    J  ~        ^=  0  J      J        ^  J 


where X (k.r)  satisfies  the  radial  equation 


(2.5) 


.2  /,_2 


(2.6) 


(dVdr'^  -  J?.  U*   l)/r  -  2M^  V°(r)/fi'^  +  k^)  X(k.r)  -  0 


%   (k.r)  ->  cos(k.r  +6.   )  as  r  ->  oo. 


Further,  the  total  energy  of  t  he  system  E  is  given  by  the  expreesion 

E  -  e,  +  ii^  k^/  2M  . 
For  our  collinear  problem  the  p-motion  is  one -dimensional,  and  the 


(H^(p)  -  ej)  JZf^(p)  ■=  0. 


vibrational  wave  fvinction  ^.ip)   is  a  solution  of  the  appropriate  vibrational 
wave  equation  - 

(2.7) 

The  cross  section  for  the  excitation  of  a  molecvile  from  a  discrete 
(bound)  vibrational  state  b  into  a  continuum  (free)  state  f  is  then  given  by 
the  following  expression  - 

(2.8) 

where 

(2.9)  V^(^)  -  I  |dp  dr  |?^(P)  ?^(k^r)  v'(p,r)  ^^(p  )  X^ik^r) . 

The  reason  for  defining  the  dissociation  cross  section  per  unit  energy  range 


d  0-^(e)/d6f  -  "  I  V^^-^  )l^  /  ^  ^f  ' 
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is  that  the  continuum  'vibrational'  wave  function  j2f„(p)  is  normalized  per 
unit  energy  range  (cf.  Section  2,d).  The  actual  cross  section  fca:  dissociation 
is  simply  obtained  by  integrating  the  expression  (2.8)  - 

(2.10)  <r^(i  ;ep  =  J  ^de^  (a  0-^^(1  jep/d  ep. 

The  character  of  the  dissociation  or  recombination  transition,  in 

which  we  pass  from  a  bound  vibrational  state  of  discrete  energy  e,  (<  O)  into 

2  2 
a  free  vibrational  state  of  energy  e„  •=  h  p„  /  2M.  (thus  6„  >  0,  i.e.,  the 

state  f  is  a  continuum  state),  should  be  examined  a  little  further.  The 

condition  of  conseinration  of  energy  is 

(2.11)  E  -  (fi^AM^)  ]^  -    \    &^    \   "  {^^  /  2M^)  k^  +  (fiV  2M^)  p^  . 

In  an  actual  situation  we  have  a  definite  boiind  vibrational  state  b  with  a 
definite  and  fixed  value  of  e,  .  The  value  of  6-  is  not  fixed,  but  covers 
some  kind  of  statistical  range  (typically  a  Maxwell-Boltzmann  distribution 
corresponding  to  the  gas  temperature).  The  incident  'external'  raomentvnn 
k-  also  covers  an  independent  statistical  range  -  which  will  of  course  also 
be  a  Maxwell-Boltzmann  distribution  for  the  appropriate  temperature  in  a 
typical  situation.  Thus,  a  characteristic  way  to  specify  energy  conservation 
is  as  follows  -  t  is  fixed;  6„  and  k„  are  treated  as  independent,  and  Eq.  (2,11) 
then  fixes  the  value  of  k,  for  a  given  set  of  values  of  6-  and  k».  To  obtain 
a  temperature-averaged  dissociation  cross  section  or  recombination  coefficient, 
we  thus  have  to  carry  out  separate  and  independent  Boltzmann  averages  over 
6»  and  k„. 


o  • 


b.  The  Translational  Potential  V  and  the  Perturbation  Potential  V 
•^^^— — — — ^-^— — — — — — —  e  -^— — ^— — — — — — — — — — — 

As  in  our  previous  work'-  *  -■,  we  neglect  the  angular  dependence  of 
V(r),  and  use  a  Morse  potential  to  describe  the  r-dependence  of  V  (r)  and  of 
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V  (p,r).  The  precise  numerical  values  are  derived  from  our  previous  v;ork  on 

m. 

V°(r)  =  U°  [e-2^(^-d>  -  e"^^^-^^]   , 


oxygen ' 


(2.12) 


< 


(2.13) 


U°  =  l.k3   X  10"-^  e^/a  :  a  =  1.1  a"-^  J  d  =  6.63  a  , 
e  'o'         0  o' 

•/   N   /      >.  ,  '   r-2a(r-z)    -a(r-z)l 

If  (p,r)  =  (p  -  p^)  U  a  |e      '  -  e      J 

u'  =  0.262  X  10"-^  e^/a  :  z  =  7.I4O  a  . 


In  deriving  Eq.  (2.13) 


'^-',  we  have  expanded  V  in  powers  of  (p  -  p  ) 


(where  p  is  the  equilibrium  nuclear  separation  of  the  AB-molecule) ,  and  have 
discarded  all  but  the  first  term.  This  has  some  justification  for  the  vibrational 
relaxation  time  problem,  since  the  m«an  amplitude  of  oscillation,  and  the  an- 
harmonicity,  of  the  lowest  vibrational  states  v  =  0,1  is  not  large.  The  expansion 
is  certainly  not  justified  in  the  present  case  on  any  quantitative  basis,  any 
more  than  is  neglect  of  the  angular  dependence  of  V(r) .  However,  it  is  certainly 
not  clear  that  putting  in  added  complexity  is  justified  by  our  knowledge  of 
intermoleciilar  potentials,  and  in  particular  these  simplifications  are  justified 
by  the  genersil  spirit  of  the  present  calculation,  which  is  to  attempt  to  get  a 
feel  for  the  general  characteristics  of  the  problem  with  a  minimal  amount  of 
numerical  work. 

c«  The  Vibrational  Motion 

3   + 
The  vibrational  levels  of  the  A  ^   state  of  0-  are  known  with 

great  accuracy  (  see  Table  I  ^  -• ) .  ITius  one  could  fit  a  Morse  potential  to 

the  dataj  however,  this  leads  to  a  rather  complicated  expression  for  the  matrix 

element  of  (  p  -  p  )  between  the  discrete  and  continuum  states.  In  view  of 

our  very  crude  approximation  to  the  perturbation  potential,  there  is  clearly 


-  8  - 


not  much  point  in  treating  the  internal  motion  in  any  very  refined  way. 
Consequently  we  have  approximated  the  potential  function  V.(p)  of  the  vibrational 
motion  by  a  square  well  of  depth  D,  the  dissociation  energy  of  the  state  A  ^  , 
and  width  L:  the  width  L  is  chosen  to  meet  the  matching  requirements  of  the 
WKB  approximation  to  the  wave  function  0,  (p)*-  -'. 

Define  ^^'^ 
(2.U4)         i?F^   /  2M^  =  D  +  e^  J  fiV  /  2M^  =  D  -  |e^|  . 

Then  for  a  square  well  of  depth  D  and  width  L,  which  is  strongly  repulsive  for 
p  <  p^    ,  the  WKB  wave  functions  are 

i^^(p)  =  N^  cos    [qCp-Pjl)  -   nA  Ji  N^  =  2A  +  1/3  ^   2/L 


(2.15) 


52f^(p)   =  N^  cos   [p(p-P;l^   ■   "/^    ]j  N^  =  2M^  /  n  fiS  , 


where  ^„(p)    is   thus  normalized  per  unit  energy  range.     For  there  to  exist  a 
boimd  state  b  of  energy  £,  ,  we  must   have   the  quantum  condition 

(2,16)  QL  =  (v  +  l/^)n  ,  V  an  integer  (vibrational  quantum  number). 


If  we  take  p,   —    p     ,  where  p     is  the  equilibrium  nuclear  separation  of  the 
AE-raolecule ,   then  the  matrix  eleme 
states  is  given  by  the  expression 


AE-raolecule,   then  the  matrix  element  of  (p  -  p  )  between  the  initial  and  final 


(2.17)  (p  -  P^)^  -  1/2  N^N^ 


(-1)^  sin  PL  _  (-l)^L  cos  PL    _    (p_q)-2 


(P-Q)^  (P-Q) 


—1  —2 

where  terms  of  order  (P+Q)       and  (P+Q)       have  been  neglected. 

d.  Flux  Conservation 

Some  of  the  results  obtained  by  the  numerical  perturbation  method 
violate  the  condition  of  flux  conservation.  That  is,  if  we  write  the  cross 
section  for  dissociation  as 
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OD 

(2.18)  <5Vk=E  (2/  +  1)      C-^ii)    , 
then 

(2.19)  (r^(i  )  >  n  /  k^^  k^  for  ^  <  ^^^.^  , 

and  n  /  k,  k     is   the  maximum  value  permitted  for  an  inelastic  scattering  cross 
section  by  flux  conservation'-    -* . 

The  fact  that  the  perturbation  method  gives  results  that  violate  flux 
conservation  is   due  to  the  relatively  small  energy  transfer  between  trans lationaL 
and  vibrational  degrees  of  freedom,  and  the  consequent  strong  coupling  between 
these  two  modes  of  motion.     For  instance,   similar  difficulties  involving 
violation  of  flux  conservation  are  qvite  familiar  in  the  case  of  elastic 
scattering  of  a  particle  by  a  center  of  force,  where  the  exact  solution  is 
of  the  form 

(2.20a)  cr(i)     =  (l^nA^)   sin^  5^      , 

while  perturbation  theory  gives  the  result^  -* 

(2.2Clb)        <r(je)  «  (UnA^)  tan^  «£  . 

A  simple  way  to  get  a  resvilt  that  may  be  meaningful  if  the  numerical 
(perturbation)  calculations  violate  flux  conservation  is  to  put 

n  /  k^^k^  for  ^  <  X  ^^.^ 


(2.21)     (r^^^(;2)  -  <; 


crP^^*  (X)  ior  Ji>Jl^^^^ 
fb 


where  J?    is  obtained  from  the  numerical  (perturbation)  cross  section 
crit 

<r  P^*  iJi)   by  using  Eq.  (2.19). 

Another  scheme  that  will  give  results  very  similar  to  those  of 
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Eq.  (2.21)  is  to  replace  it  by 


(2.21*)   (T^^^  {£)  i 


0        for  i  >^ 


crit 


e.  Three-Body  Recombination  Coefficient 

When  two  atoms  A,B  combine  from  the   'free'   state  f  of  the  molecule 
AB  into  the  boiind  state  b,  in  the  presence  of  a  third  body  C   (to  conserve 
energy  and  raomentura),  the  nvunber  of  recombinations  per  unit  volume  and  iinit 
time  is  given  by 

^2-22)  d  n^^^  /dt  =  C^f  n^  n^  n^,     , 

where  n.,  Ur,,  n-  are  the  munbers  of  A,B,C  atoms  per  unit  volume;  C,  -  is  the 
three-body  recombination  coefficient. 

Three-body  recombination  is  the  inverse  process  to  dissociation  of 
the  molecule  AB  by  collision  with  C,  and  thus  the  three-bocfy  recombination 
coefficient  C,  „  is  related  to  the  dissociation  cross  section  C"^  by  the 
principle  of  detailed  balance*  In  particular,  for  the  reaction  (1.2),  the 
number  of  dissociations  per  unit  volume  per  unit  time  is  given  by 

(2.23)        d  n^33  /  dt  =  o-^  {i\MJ   n^,  n^„    . 

In  a  state  of  statistical  equilibrium, 

(2.2U)         d  n^33  /  dt  .  d  n^^J  dt  , 

so  that 

Sf  "  (^AB-  ^C"/  ^A  "b  V  (*^/"e^  ^fb   • 

ri7"i 

Now,  from  equilibriiam  considerations'-  -*   we  obtain  the  result 


-  n  - 


(2.26) 


2  2 


(^AB-  "C"  /  "a  "b  "c^  "  ^M  '   ^AB^A  °B^^*^AB^  ^^"  %    ^/^fV  ' 


which  is  derived  in  Appendix  A,  For  ovir  process  (l.l)  we  have  6.  =  1,  and 


AB 


G 


^  =  3  G^  Gg  ,  so  that 


(2.27) 


Tbf  »  6  n2  (  p2  k^  /  p^  icp 


We  are  interested  in  the  temperature-averaged  three-body  recombination 
coefficient.  From  the  disciission  below  Eq.  (2.11)  we  see  that  before  the 
temperature  averaging,  C,  _  depends  on  any  two  of  (k,  ,  k„,e„),  because  of  the  energy 
conservation  relation  linking  these  quantities.  For  definiteness,  let  us  treat 
k-  and  6-  as  independent  variables.  Under  normal  conditions,  k„  and  6-  are 
separately  and  independently  distributed  in  Maxrrfellian  fashion  characterized 
by  a  temperature  T.  Thus,  integrating  first  over  e„,  we  get 

1/2  ^-6„/kT 


'  ^f  (  kf^  ^f,avT>  'l^'f    Tbf  («S/^e^  ^fb  ^^f'^f^  ^  'Y^  ''"'' 


(2.28) 


A 


2n  /  (nkT) 


3/2 


and  finally,  carrying  out  the  second  Maxwell-Boltzmann  average,  we  get 

.2,2 


'   C^^(T)   -  A*     J     k^  dk^  e-^  ^f  /  2^e  ^     C^,  (k,,   e,^^^) 


(2.29) 


A* 


Ijn   (fi^  /  2nM  kT  )^/^     . 


The  integrations  extend  over  0<6.  <oo,0<k^<a);  o^  Jvist  vanishes  for 
those  valxies  of  (e„,  k.)  that  do  not  satisfy  energy  conservation. 

The  actual  ntimerical  calculations  were  nade  by  averaging  over  e^  and 
k.  -  not  k-  -  with  independent  Maxwellian  distributions.  In  practice,  if  one 
takes  energy  conservation  into  account  correctly,  this  makes  no  essential 
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difference,  and  the  numerical  results  in  the  present  situation  are  certainly 
not  affected.  However,  on  physical  grounds  it  seems  better  to  average  over 
k^  and  &     rather  than  over  k,  and  e-  . 

f •  A  Simple  Physical  Model  of  the  Recombination  Coefficient 

In  a  situation  where  perturbation  theory  fails  because  the  predicted 

cross  section  is  so  large  that  it  violates  flux  conservation,  we  can  try  to 

get  a  rough  analytical  estimte  of  ^   ..  of  Eq.  (2,19),  and  then  use  Eq.  (2, 21) 

to  give  a  value  of  ov_  and  thus  get  C,  „. 

Evidently  if  the  intermolecular  potential  has  some  mean  range  R,  then 

this  gives  us  a  critical  angular  momentum  j?   ..  from  the  standard  'impact 

parameter'  condition  - 

(2.30)         *^crit"*^^  • 

If  we  make  use  of  the  fact  that  k„  —  k.    and  substitute  this  and   (2.30)  in 
(2.21),  we  get  the  result 

X"   0 
and  substituting  this  in  Eq.   (2.28),    (2.29),  we  get 

(2.32)  ^bf^""*^^   '  ^^"^  Po  «^  ^  /  (Mi^g)-^/^     . 

We  notice  that  the  result  (2.32)  as  written  down  is  a  fvinction  of 

temperature  T,  only  insofar  as  it  affects  the  effective  range  R,  and  of  the 

vibrational  state  v  only  to  the  extent  that  it  enters  in  p  *  p  (v).  We  do 

■^o    o 

observe,  however,  that  this -admittedly  very  rough-approximation  indicates 
that  C,  -  increases  with  increasing  v  (which  increases  the  effective  p„(v)  ), 
and  decreases  with  increasing  temperature  since  the  atoms  can  approach  closer 


We  should  like  to  thank  Dr.  K.  Takayanagi  for  pointing  this  out  to  usf 
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and  closer. 

To  get  a  simple  ntunerical  estimate,  let  \xs  put  R  =  7  a     —   d  of 
Section  2.b,  then  from  Table  I  we  have  p  (v)  for  v  -  9,10,  which  gives 


(2.33) 


C^^^PP'(v  -  9)   -  6.I4  X  10"^^  craVsec 


C^^^PP*(v«10)   -  8.0  X  10"^^  cmVsec   . 

A  number  of  approximate  estimates  of  three-body  recombination 
coefficients  have  been  made  by  different  people.     The  simplest  sort  of  con- 
siderations  -  obtained  virtually  from  dimensional  analysis  -  give*-    '       -* 

(2.3U)  C^gj,  -  f  R^  V     , 

where  v  =  (SkT/nM)''''^'^  is  the  mean  speed,  M  =  reduced  mass,  and  f  is  some 
kind  of  a  probability  factor.  For  f-^l,  R  =  7a,v»  10  cm/sec,  this  gives 

(2.35)         C   -~  5  X  10"-^^  cmVsec  . 
rec 


3.  Results 

The  results  of  the  Univac  calculation  are  given  in  Tables  2  and  3. 

ri9  20I 
They  are  compared  in  Figure  1  with  the  latest  experimental  results  '-*-'. 

We  see  that  the  experimental  and  theoretical  results  agree  to  within  a  factor 

of  3,  which  is  a  very  considerable  improvement  over  previous  detailed  quantum 

mechanical  calculations  L  >  MJ .  ctae  thing  that  is  demonstrated  quite  clearly 

by  the  calculations  is  that  recombination  will  lead  to  vLbrationally  excited 

molecules.  This  has  been  predicted  by  the  theorists  for  qiiite  a  long  time"-  *  ^J, 

but  thus  far  there  is  not  very  much  definite  experimental  evidence  for  it*-  -'. 

To  be  more  specific,  it  has  alMays  been  hard  on  strictly  theoretical 


grovinds  to  understand  the  large  experimental  values  cf  recombination  rates j 
the  fact  that  the  present  calculations  -  which  take  due  account  of  the 
limitations  on  cross  sections  imposed  by  flux  conservation  -  still  exceed 
the  experimental  values  is  thvis  an  error  in  the  right  direction.  This  may- 
be corrected  by  a  more  appropriate  choice  of  intermolecular  potential,  perhaps, 
or  by  detailed  considerations  that  might  indicate  that  a  certain  fraction  of 

0  (  P)  atoms  will  recombine  into  some  state  such  as  Op  (X  ^  ),  which  might 

g 
perhaps  have  a  state  with  binding  energy  significantly  greater  than  the  state 

Op  (a  ^    ;  v  =  10),  thus  giving  a  smaller  effective  recombination  coefficient. 

u 
Of  course  it  must  be  stressed  that  the  experimental  data  are  subject  to  some 

errors,  while  the  theoretical  results  certainly  cannot  be  relied  upon  to  within 
better  than  a  factor  of  order  3,  on  accovmt  of  -uncertainties  arising  from  such 
considerations  as  the  intermolecular  potentials  or  detailed  knowledge  of  boiind 
vibrational  states.  A  further  source  of  error  in  the  present  calculation  is 
that  the  effects  of  molecular  rotation  are  neglected  and  all  collisions  are 
treated  as  though  the  potential  were  the  same  as  for  a  collinear  collision  - 
an  obvious  overestimate. 

To  sum  up,  it  was  clear  fran  the  outset  that  a  calculation  of  a  problem 
as  complicated  as  three-body  recombination  could  not  aim  to  obtain  detailed 
quantitative  results  because  of  all  the  various  factors  that  enter  when  one 
really  tries  to  take  all  relevant  physical  effects  into  accountj  and  to  make 
an  estimate  of  these  -  even  if  one  were  willing  to  undertake  the  labor  -  would 
call  for  detailed  experimental  information  that  is  simply  not  available.  The 
present  calculation  is  thus  a  schematic  one  using  a  greatly  oversimplified 
model.  However,  within  its  obvious  limitations  the  present  calculation  does 
indicate  that  the  orders  of  magnitude  obtained  in  the  experiments  can  be  ex- 
plained in  terms  of  'conventionally  reasonable'  mechanisms. 
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There  is  one  additional  piece  of  information  that  can  be  obtained, 
and  that  relates  to  the  temperature  dependence  of  the  recombination  coefficient. 

ri9i 

We  see  that  the  experimental  results'-  -"  show  very  little  temperature  dependence. 
The  calculated  recombination  rate  into  the  state  v  =  10  does  not  vary  veiy 
strongly  with  temperature,  but  still  there  is  a  low  temperature  fall  in  the 
rate.  A  more  uniform  rate  would  result  if  one  supposed  that  recombination 
occurs  into  a  state  with  binding  energy  lower  than  the  100  cm   that  has  been 
used  here  for  the  state  v  =  10,  or  perhaps  into  several  such  states.  It  might 
just  be  possible  to  obtain  information  on  such  very  weakly  bound  vibrational 
states  from  some  considerations  related  to  these. 
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Appendix  A.  Detailed  Balance 

For  the  purpose  of  the  present  discussion,  which  seeks  to  obtain 
a  relation  between  the  dissociation  cross  section  <J^  and  the  three-body 
recombination  coefficient  C,  „,  we  consider  that  the  raoleciole  AB  is  formed 
with  rotational  energy 

(A.i)       E^^^  =  C-h^/aM^  Q^   j(j  +1)  . 

This  is  very  convenient  for  our  present  purpose,  and  does  not  conflict  sig- 
nificantly with  the  spirit  of  the  perturbation  calculation,  even  though  we 
only  consider  collinear  collisions  there.  Thus,  we  may  define  E.  as  the 
net  energy  transferred  between  the  'internal'  and  'external'  degrees  of 
freedom,  and  then  conservation  of  energy  gives  the  two  relations 


(A.2) 


(a)    E,^  =  (fiV2M^)  (1^  -  k^) 


(b)    Eie  =  (*V2M.)p/-E^^^*  1.,!  . 

Now,  both  the  dissociation  cross  section  0%.,  and  the  recombination 
coefficient  C,  »  are  purely  atomic  properties  which  are  independent  of  any 
macroscopic  properties  of  the  system,  such  as  whether  or  not  it  is  in  equi- 
librium, of  the  temperature,  density,  etc.j  but  just  because  of  this  fact 
we  may  obtain  the  detailed  balance  relation  under  equilibrium  conditions, 
which  makes  the  calculation  quite  simple*-  -■.  While  the  temperatvire  T,  and 

explicit  references  to  equilibrium  conditions,  enter  into  the  course  of  the 
analysis,  these  quantities  fall  out  before  the  final  result. 

ri7i 

Equilibrium  considerations  give  the  result '-  -^ 
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f  ^Pf 


v;here  n.     ,  n^  °  ,  etc.  refer  to  the  total  number  densities  of  A,E,  etc.j 
n,,  Hp  etc.  refer  to  the  number  densities  for  momenta  p„,  k.,  etc.  f,(T,M)  is 
the  translational  partition  function  for  a  syston  of  mass  M  at  temperature  T, 
M  =  (M.+  Vi^)   f  (T)  is  a  rotational  and  f-(T)  a  vibrational  partition  function. 
The  term  (2J  +  l)  arises  from  the  rotational  state  of  the  molecule;  h  is  the 
balance  in  ths  phase  space  volvime  factors  arising  between  (AB  +  C)  on  the  one 
hand,  and  (A+E+C)  On  the  other.  The  other  terms,  (Ijnh  k„  dk^),  etc.  are  phase 
volume  contributions. 

We  have  assumed  no  change  in  the  electronic  state  of  the  system.  The 
total  n\mber  densities  n,  '°  ....  are  related  in  equilibrium  by  the  reaction 
isochore  <-  J  _ 

(tot) 

i 

^A     "b 


°AB  ^1^^'^AB^  ^2^"^^  ^3^'^ 


/.  ,s        "ab  

^  ^^       n/tot)   (tot)   ^U-fi^B^'^A^B    •   V^^^WV^A^ 


where  G. ,  G_,  G._  are  the  electronic  statistical  weights  of  A,B,  AB,  and 

6   =  1  if  A  and  B  are  identical,  and  zero  otherwise  (nuclear  symmetry  factor). 

Substituting  Eq.  (A.I4)  in  (A. 3)  and  cancelling  out  tenns,  we  get  the 

result  - 

2 


^-  2n^   (2J.^  \^ 

bf    G.  G„  (1  +6.^;   '     *   2 


To  get  rid  of  the  differentials  that  arise  from  the  phase  space  volume,  we 
differentiate  Ec^.  (A. 2a),  (A. 2b),  keeping  E.   and  6^  constant;  thus  Eq.  (A. 2a) 
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and  (A. 2b)  give,  respectively, 
(A.6a)        kj^  dkj^  =  k^  dk^  , 
and 
(A.6b)         p^  c^^  -  (2J  +  1)  2  p^  . 

Notice  that  to  get  two  relations  (A. 6a),  (A. 6b)  from  the  single 
condition  of  conservation  of  energy  (2.11)  we  have  made  an  additional  assumption 
implicit  in  Eq.  (A. 2)  and  explicit  in  the  differentiation.  That  is,  the  phase 
space  surface  over  which  the  condition  of  detailed  balance  is  applied  is  chosen 
so  that  E.  is  constant  along  it. 

If  one  now  substitutes  Eq.  (A. 6)  in  (A. 5),  one  gets  the  final  result 


-  IS   - 


Appendix  B»  Dissociation  Rates 

From  our  aralysis  it  is  quite  easy  to  estimate  the  temperature- 
averaged  dissociation  rate  k,.      (T),  which  is   just  the  temperature  average 
of  (file /M  )   <T„»     The  total  rate  of  dissociation  is 

(B.l)         (^M)  „,,^^  -  k^^^(T)  n^'*°*'  n,(*°*)       . 


ri7~i 

From  Eq.   (2.23),  and  the  standard  result"-    -■ 


2,  2 


/.,   ^N  /,  N    .,  (tot)      (tot)   ,^,2    ,,        fi^k,    /2M  kT-(D  -k  DAT 

(B.2)         n^,n^„   (k^,^)   dk^  =n^^        ^n^  ^*  ^S  ^'S  ^       ^         «  °^         > 

where  the  constant  A*  comes  from  Eq.  (2.29)  and  D  is  the  dissociation  energy 
of  the  ground  state,  we  get  the  result 

Notice  that  we  have  used  <r^  (k^,  e^  avT^»  which  is  given  by  the 
expression  - 

(E.M    cr^  i)^,   6^  ^^)  =  A     cr^  (k^,  e^)  e/  e  f^  de^, 

^0 


where  A  is  defined  in  Eq.  (2,28). 

It  is  quite  easy  to  get  an  expression  for  k^j^gg(^)  i"  terms  of  the 

approximation  scheme  of  Section  2.f,  which  already  gives  us  effectively  a  value 

of  <r.  (k^,  6„  J   -  cf.  Eq.  (2.31)  for  instance.  Thus  we  get  the  result 

(B.5)    k^ss^^^*^^^  "  "  R^(8kT/TiM^)^/2  (1  +  |e^|/kT)  e'^o^   . 

Note  that  this  result  has  been  worked  out  for  the  case  when  we  carry  out  independent 
Boltzmann  averages  over  e^  and  k^^,  rather  than  over  e^  and  k^.   In  practice 
this  should  not  make  any  significant  numerical  difference  to  the  results* 
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3  ^  + 


Table  1.   Highly   \abrationally  Excited  Levels  of  the  A"^  ^^ State  of  0^ 

vibrational  quantum  v  binding  energy   |6,(v)|(cm~  )  nuclear  separation  Pp(v)(a^) 
10                                                      lOli  3.69 

9  3liO  3.56 

••••  ••*•  •••• 

0  5570  2.90 
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Table  2«  numerical  Results  for  the  Capture  Cross  Section  <r^(e.) 
Vibrational  state  v  ■  10  Vibrational  state  v  =  9 


h 

^f 

-^crit 

^fb^^f^ 

15 

1 

— 

26 

2 

ao 

81 

3 

60 

93 

U 

99 

5 

> 

f 

105 

20 

1 

— 

53 

2 

80 

111 

3 

(>100) 

121 

a 

12U 

5 

127 

6 

> 

f 

132 

25 

1 

60 

72 

2 

120 

112 

3 

(>ll40) 

135 

u 

137 

5 

> 

'• 

139 

30 

1 

100 

73 

2 

iho 

89 

3 

(>160) 

93 

ii 

9U 

5 

> 

r 

95 

\ 

''f 

0 
'^crit 

^fb^^P 

20 

1 

— 

2.I4 

2 

— 

8.1 

3 

— 

12. a 

U 

— 

15.2 

25 

1 

— 

7.6 

2 

— 

27 

3 

Uo 

ii2 

a 

60 

51 

30 

1 

15.5 

2 

80 

a6 

3 

100 

52 

U 

100 

5a 

5 

(>120) 

56 

35 

1 

— 

32 

2 

120 

68 

3 

150 

80 

l4 

150 

81 

5 

150 

82 
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Table   2  -  Continiied 


V  =  10 

■s 

"■t 

-^crit 

'^fb^^P 

liO 

1 

180 

8U 

3 

210 

108 

5 

210 

117 

7 

(>2l40) 
r 

118 

9 

121 

V  =  9 


^ 

^f 

-^crit 

•^fb^^p 

llO 

1 

120 

^2 

3 

180 

86 

5 

93 

7 

Y 

96 

50 

1 

180 

68 

3 

2ijO 

100 

5 
7 

(>270) 

117 
119 

Notes 

—1  2 

(i)    the  units  of  k,  ,  P-  are  a   j  cr™  (c.)  is  measured  in  a  . 

(ii)   for  the  definition  of  ^  ..  see  Eq.  (2.19);  see  Eq.  (3.1)  for  <^^(^f) 

when  there  exists  a  value  of  ^  ...  i.e.  when  flux  conservation  limits  the 

crit 

cross  section. 

(iii)  a  numerical  value  of  ^  ..  in  parentheses,  e.g.  (>  100)  means  that  flux 

crio 

conservation  was  violated  up  to  the  highest  value  of ^  for  which  calculations 
were  made.  Effective  values  of  <r„  (e„)  could,  however,  be  inferred  from  other 


results. 
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Table  3.   Calculated  Values  of  Three-Body  Kecombination  Coefficients  C^^(T;v) 

,^  /      N  ^     fn.  \        I     ^  /       )  C^„(T;   V  -  9)      (cm  /sec) 

kT  (ev)  ^bf^    '   ^  "  ^  ""  ^^®° 


03 

(2.5  -  3.5)10"^^ 

10 

5.6  X  10"^^ 

20 

6.3  X  10"^^ 

30 

6.3  X  10"^^ 

1.1  X  10"^^ 
2.3  X  10"^^ 

3.2  X  10"^-^ 


Note  -  the  values  at  kT  =  .03  ev  (350°K)  are  uncertain  because  there  were 
not  enough  numerical  results. 
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